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Ultra-cold bosons in zig-zag optical lattices present a rich physics due to the interplay between 
frustration, induced by lattice geometry, two-body interaction and three-body constraint. Uncon- 
strained bosons may develop chiral superfluidity and a Mott-insulator even at vanishingly small 
interactions. Bosons with a three-body constraint allow for a Haldane-insulator phase in non-polar 
gases, as well as pair-superfluidity and density wave phases for attractive interactions. These phases 
may be created and detected within the current state of the art techniques. 



Introduction Atoms in optical lattices offer extraordi- 
nary possibilities for the controlled emulation and analy- 
sis of lattice models and quantum magnetism [l[ . Various 
lattice geometries are attainable by means of proper laser 
arrangements, including triangular Q and Kagome Q 
lattices, opening fascinating possibilities for the study of 
geometric frustration, which may result in flat bands in 
which the constrained mobility may largely enhance the 
role of interactions Moreover, the value and sign of 
inter-site hopping may be modified by means of shak- 
ing techniques , allowing for the study of frustrated 
antiferromagnets with bosonic lattice gases 0. 

Interatomic interactions may be controlled basically at 
will by means of Feshbach resonances In particu- 
lar, large on-site repulsion may allow for the suppres- 
sion of double occupancy in bosonic gases at low fill- 
ings (hard-core regime). Interestingly, it has been re- 
cently suggested that, due to a Zcno-like effect, large 
three-body loss rates may result in an effective three- 
body constraint, in which no more than two bosons may 
occupy a given lattice site Q. This constraint opens 
exciting novel scenarios, especially in what concerns sta- 
ble Bose gases with attractive on-site interactions, in- 
cluding color superfluids in spinor Bose gases [llj and 
pair-superfluid phases 0, [HI, 12 1. The suppression of 
three- body occupation has been hinted in recent experi- 
ments [13j . 

Under proper conditions, lattice gases may resemble 
to a large extent effective spin models, e.g. hard-core 
bosons may be mapped into a spin-1/2 XY Hciscn- 
bcrg model Lattice bosons at unit filling resemble 
to a large extent spin-1 chains [l4| . and in the pres- 
ence of inter-site interactions, as it is the case of polar 
gases (lif , have been shown to present a gapped Haldane- 
like phase [HI (dubbed Haldane-insulatoWHI) [lUEzi) 
characterized by a non-local string-order [1 81 ] . 

In this Letter we analyze the physics of ultra-cold 
bosons in zig-zag optical lattices. We show that the 
interplay of frustration and interactions lead to a dif- 
ferent physics for unconstrained and constrained (with 
up to two particles per site) bosons. For unconstrained 
bosons, geometric frustration induces chiral superfluid- 
ity, and allows for a Mott-insulator phase even at vanish- 
ingly small interactions. For constrained bosons, we show 




FIG. 1: Zig-zag chains formed by an incoherent super- 
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V 2 {r) = V 20 sm 2 (V3ky/4 - tt/4). In the fi gure, in which 
V20/V10 = 2, darker regions mean lower potential. The 
hopping rate between nearest (next-nearest) neighbors is 
t > (*')■ 

that a Haldane-insulator phase becomes possible even 
for non-polar gases. Moreover, pair-superfluid [l, [H, 12| 
and density-wave phases may occur for attractive on-site 
interactions. A direct first-order phase transition from 
Haldane-insulator to pair-superfluid is observed and ex- 
plained. These phases may be realized and detected with 
existing state of the art techniques. 

Zig-zag lattices. In the following we consider bosons 
in zig-zag optical lattices. As shown in Fig.Q] this partic- 
ular geometry may result from the incoherent superpo- 
sition of a triangular lattice with elementary cell vectors 



Si = ae x and 0*2 



e y ) (formed by three laser 



beams of wavenumber k = 47r/3a oriented at 120 degrees 
from each other, as discussed in Ref. 0]) and a superlat- 
tice with lattice spacing y/3a oriented along y. For a suf- 
ficiently strong superlattice, zig-zag ladders are formed, 
and the hopping between ladders may be neglected. We 
will hence concentrate in the following on the physics 
of bosons in a single zig-zag ladder, which is to a large 
extent given by the rates t and t' characterizing the hop- 
ping along the two directions a\^. (Fig. [IJ. As shown in 
Ref. 0, a periodic lattice shaking may be employed to 
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control the value of t and if independently. Interestingly, 
their sign may be controlled as well. In the following we 
consider an inverted sign for both hoppings, which result 
in an anti- ferromagnetic coupling between sites Q- 

Model. Ordering the sites as indicated in Fig. [TJ the 
physics of the system is given by a Bose- Hubbard Hamil- 
tonian with on-site interactions characterized by the cou- 
pling constant U, nearest-neighbor hopping t < and 
next-nearest-neighbor hopping if < 0: 
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where b\,bi are the bosonic creation/annihilation opera- 
tors of particles at site i, rii = b\bi, and we have added 
the possibility of three-body interactions, characterized 
by the coupling constant f/ 3 . We assume below an aver- 
age unit filling n = 1. 

Unconstrained bosons. We discuss first the ground- 
state properties of unconstrained bosons (C7 3 = 0). At 
[7 = 0, the Hamiltonian (TTJ) is diagonalized in quasi- 
momentum space H = e(k)b\bk, with the dispersion 
e(k) = |f|(cosfc + j cos2fe), with j = t'/t. Depending 
on the frustration j we may distinguish two regimes. If 
j < 1/4, the dispersion e(k) presents a single minimum 
at k = 7r, and hence small U will introduce a super- 
fluid (SF) phase, with quasi-condensate at k = tt. If 
j > 1/4, e(fc) presents two non-equivalent minima at 
k = kO = ± arccos[— l/4j]. As shown below, interac- 
tions favor the predominant population of one of these 
minima, and the system enters a chiral supcrfluid (CSF) 
phase with a non-zero local boson current characterized 
by a finite chirality (k^), with Hi = ^{b\bi + \ — H.c). At 
j = 1/4, the Lifshitz point, the dispersion becomes quar- 
tic at the k = it minimum, e(k) ~ {k — 7r) 4 , the effective 
mass m = (d 2 e(k) / 'dk 2 )^^ = l/t(l — 4j) diverges, and 
even vanishingly small interactions become relevant. 

To study the effect of interactions we combine numer- 
ical calculations based on the density matrix renormal- 
ization group (DMRG) method ^ (with up to N = 300 
sites keeping per block on average 400 states for gapped 
phases and 600 states for gapless ones), and bosoniza- 
tion techniques to unveil the low-energy behavior of 
model (TTJ). For j < 1 /4, we employ standard bosoniza- 
tion transformations [19| . with an additional oscillating 
factor bi — s- (— iy e l V^ e (x) ^ ^ Q obtain the low-energy effec- 
tive theory, which is given by the sine-Gordon model 
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-M cos [2imx — V 47T0] , (2) 



where 6 and d x <p describe phase and density fluctuations 
of bosons respectively, [9(x),d y (j)\ = i8(x — y), v s is the 
sound velocity and K the Luttinger parameter. In the 
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FIG. 2: Phase diagram for unconstrained bosons as a function 
of the frustration parameter j and (a) the on-site interaction 
U (with U3 = 0) and (b) the three-body repulsion U3 (and 
U = 0). In the figures, O indicate the boundary of the chiral 
phases characterized by long-range ordered chirality-chirality 
correlations {niftj); □ indicate the boundary of the SF-phases 
indicated by the critical Luttinger parameter K = 2. Note 
that narrow CMI and CHI phases may occur as well. For 
U, U3 > 0.5 keeping n max = 4 bosons per site in the DMRG- 
simulation can be shown to be sufficient. 



weak-coupling, Um 1, hydrodynamic relations are ex- 
pected to hold: v s (j) ~ \JnU '/run 2 = v s (0)- v / l — 4j and 
K(j) ~ yj niT 2 /Um = K(Q)y/l — 4j, clearly showing that 
j enhances correlations. At j = 1/4, m diverges and the 
system enters a Mott-insulator (MI) even for vanishingly 
small U (Fig. EJa)). The SF-MI transition takes place 
however in the strong-coupling regime in which v s and 
K must be determined numerically. We obtain K from 
the single-particle correlations Gij = (bibj) which in the 
SF decay as - - j\~ 1/2K ■ The value K = 2 

marks the boundary between SF (U < U c , K > 2) and 
MI (U > U c , K < 2, and M > 0). The MI phase 
is characterized by a hidden parity order [T3] , 2 P = 
lini| i _ i |^ 00 ((-l)^'< i <j 5ni ) - (cos y/ircf)} 2 , which has 
been recently measured in site- resolved experiments [2(| ■ 
The j > 1/4 case is best understood from bosonization 
in the j > 1 regime. We may then introduce two pairs 
of bosonic fields (0i,<fii) and (02,</'2), describing, respec- 
tively, the subchains of even and odd sites. The effective 
model is governed by the Hamiltonian density 
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\d x 0+ sin v 2it9- — 2M. cos v 27r</> + cos \2nij) 



where 9 ± = {9 1 ± 6 2 )/V2, <f>± = {4>i ± fa) 1^/2, v±, K ± , 
and M. are phenomenological parameters (in the regimes 
displayed on Figure 2 (a)), A ~ j -1 . Note that the chi- 
rality is given by — > sinv / 27r^-(a;). In weak-coupling, 
Um' < 1, with to' = {d 2 e{k)/dk 2 )^ = 4j/i(16j 2 - 1), 

v± ~ \/nU /m'ir 2 and K± ~ y 'nn 2 /L r m' '. In this case 
only the term d x Q+ sin \[2/kQ^ is relevant, resulting in 
(sin 1 



2tt9-) 7^ [2l|. Hence, a small U is expected 



to favor a CSF for j > 1/4, as our numerical results 
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confirm (Fig. HJa)). The CSF phase is characterized by 
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Moreover, depending on the values of K± bosoniza- 
tion opens the possibility of two consecutive phase transi- 
tions with increasing U starting from the CSF phase [22| , 
which we have confirmed with our DMRG calcula- 
tions (Fig. [Sfa)). First a KT transition occurs from CSF 
to chiral-Mott (CMI), a narrow Mott phase with finite 
chirality. Then an Ising transition is produced from CMI 
to non-chiral MI. At both KT transition lines in Fig.[2Ja) 
(SF-MI and CSF-CMI), up to a logarithmic prcfactor, 
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where in CSF « f 0. 
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Constrained bosons. As mentioned above, sufficiently 
large three-body losses may result in a three-body con- 
straint (&f) 3 = (U 3 = oo) @. In that case, Model JT]) 
may be mapped to a large extent onto a frustrated spin- 
1 chain model |23j |. which, presents the possibility of 
a gapped Haldane phase, characterized by a non-local 
string order. Hence, interestingly, constrained bosons in 
a zig-zag lattice may be expected to allow for the observa- 
tion of the HI phase in the absence of polar interactions. 

Indeed, a model with [7 = and finite U3 
shows that at the Lifshitz point, j = 1/4, a HI 
phase is stabilized for arbitrarily weak U3 (Fig. HJb)). 
The effective theory describing the HI is again the 
sine-Gordon model © with K < 2. However, 
now J\4 < 0, which selects a hidden string or- 
der 0% = ]im.\ i -j\- yoo (6mexp[i'X J2i<i<j ■ SnijSnj) ~ 
(sin 0F(9) 2 Resembling the case of Fig.^a), SF, HI, 
chiral-HI (CHI) and CSF phases occur (Fig.Hfb)). Th ese 
phases are expected for U3 = 00 from known results in 



frustrated spin- 1 chains [27H29|. Our DMRG simulations 
suggest that all these phases meet at j = 1/4 for U3 —> 0. 

Figure |3] shows the phase diagram for constrained 
bosons (U3 = 00). Starting from the HI phase, increas- 
ing U > can induce a Gaussian HI-MI phase transition, 
characterized by a vanishing J\A = in ([2]), resembling 
the phase transition between Haldane and large-D phases 
induced by single- ion anisotropy in spin-1 chains (30| . 
The SF phase is separated from the MI and HI by 
KT transitions, whereas at the CSF boundary with the 
MI (HI) a CMI (CHI) occurs as mentioned above (these 
very narrow regions are not resolved in Fig. [3]). 

Interestingly, constrained bosons allow as well for the 
exploration of attractive two-body interactions, U < 0, 
without collapse. The U < phases are also depicted in 
Fig. El For sufficiently large \U\, bosons tend to cluster 
in pairs and, as already discussed in Ref. |{|, for j = an 
Ising transition between a SF and a pair superfluid (PSF) 



occurs 



24j, analogous to the XY1 to XY2 phase transi- 
tion in spin-1 chains induced by single- ion anisotropy [30j ] 
(this transition has been recently studied for 2D lattices 
as well [III EH)- The PSF phase is characterized by 
an exponentially decaying Gij but algebraically decay- 
ing pair-correlation function g\ 2 = ((b\) 2 (bj) 2 ) ■ Indeed 
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FIG. 3: Phase diagram for constrained bosons as a function 
of the on-site interaction U and the frustration parameter j. 
Narrow CMI and CHI phases may occur along the phase tran- 
sition lines from CSF to MI and CSF to HI, respectively, but 
their extension would be negligible in the figure. For the pre- 
cise location of the PSF-DW and HI-MI transition we have 
additionally analyzed the energy-level-crossings with, respec- 
tively, periodic and twisted-boundary conditions. The KT 
transitions from SF to MI and HI have been determined by 
the extraction of the Luttinger parameter K [26l |. 

a PSF occurs for sufficiently large \U\, also for j > 1/4 
which is characterized in bosonization in Eq. ([3]) by a 
gapped antisymmetric sector, with pinned </>_, and a gap- 
less symmetric sector |3lj |. Though one may anticipate 
an Ising phase transition between the CSF (with broken 
discrete parity symmetry) and the PSF (with restored 
symmetry), the behavior of 2 P and k (not shown) hints 
to a weakly first-order nature. 

Small U < disfavors singly-occupied sites and 
thus enhances C| and the bulk excitation gap of the 
HI phase (see Figs. 2] and [S]). However, since large 
U < removes singly occupied sites completely, just 
like strong nearest neighbour repulsion, it is expected 
that the HI phase eventually will transform for growing 
\U\ into a gapped density- wave (DW) phase via Ising 
phase transition [lij], and string order will evolve into 
DW order (Fig. 2] shows how O s merges with Odw = 



lira., 



>(-l) J ( 



TliUi+j, 



for U/t < -3). The DW phase is 



characterized by an exponential decay of both Gij and 

(2) 

G\j' though a finite Odw- Our DMRG results confirm 
this scenario (see Fig. 2]) , showing that a DW phase is lo- 
cated between the above mentioned PSF regions (Fig. [3]) . 

Interestingly the DW phase remains in between both 
PSF regions all the way into U — s- — 00. In that regime, 
we may project out singly-occupied sites, and introduce 
a pseudo-spin- 1/2, identifying |0) — > \]), |2) — > |f), and 
defining the spin operators rr — > (— iyb 2 /^/2, 2rf — > 
b\bi — 1. The effective model to leading order in 1/\U\ is 
a spin-1/2 chain: 
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FIG. 4: Order parameters 2 P (□), 0% (A), O dw (v), and 
energy gap (x) as function of U for constrained bosons on 
the j — 0.3 line (N — 160). The parity, as defined in the 
Mott state, must get an additional oscillating factor in the 
DW phase 2 P (-l) 4 " J Op. 



FIG. 5: Order parameters Op (□), 0% (A), and k 2 (O) as 
function of U for constrained bosons on the j = 0.6 line (N = 
160). Both 2 P}S show jumps at the HI to PSF transition 
at U ~ -2.5. In the PSF, similar to the DW phase, 2 P is 
defined with an additional oscillating factor. 



where J = t 2 /\U\. For j = 0, this is a SU(2) symmet- 
ric chain, whereas the j 2 J terms break the symmetry 
down to U(l), moving the effective theory obtained after 
bosonization of Hi towards the irrelevant direction (in 

2 v 

the renormalization group sense). As a result of this, a 
gapless XY phase of the spin-1/2 chain is expected, i.e. 
a PSF phase. Higher order terms in 1/\U\ (not shown 
explicitely) break, even for j = 0, the SU (2) symmetry 
to U(l) in the irrelevant direction. However, interest- 
ingly, the ring exchange along the elementary triangle of 
the zig-zag chain, with amplitude jt 3 /U 2 , forces the ef- 
fective theory towards the relevant direction, leading to 
a gapped Neel phase of the spin-1/2 chain, i.e. the DW 
phase. The competition between exchange along the lat- 
tice bonds and ring-exchange leads hence to two consec- 
utive KT phase transitions induced by j, for j<l first 
from PSF to DW, followed by DW back to PSF. The 
width of the DW phase is ~ t/\U\, and it extends all the 
way into the U — > — oo limit. 

Finally, our DMRG simulations show a narrow re- 
gion where a direct, apparently first-order, HI-PSF tran- 
sition occurs, characterized by discontinuous jumps of 
Og P (Fig. [5]). This first-order nature is explained be- 
cause on one hand increasing \U\ within the HI phase 
increases C| due to the suppression of singly-occupied 
sites, and on the other hand, for 0.6 < j < 0.75 (Fig. [3]), 
a growing \U\ destroys the insulating state in favour of a 
PSF phase, where string order cannot exist. On the con- 
trary, Og diminishes for decreasing \U\ when approaching 
the HI-CSF boundary (Fig. E) 

Conclusions. In summary, the interplay between geo- 
metrical frustration and interactions leads to rich physics 
for ultra-cold bosons in zig-zag optical lattices. Un- 
constrained bosons may present chiral superfluidity, and 
Mott insulator for vanishingly small interactions. Con- 
strained bosons may allow for the observation of Haldane- 
insulator without the necessity of polar interactions, as 
well as pair-superfiuid and density wave phases at at- 
tractive interactions. All the predicted phases may be 
detected using state of the art techniques. The SF and 
CSF phases may be distinguished by means of time-of- 
fliglit (TOF) techniques, in a similar way as recently done 



for condensates in triangular optical lattices 0| • The DW 
and PSF phases are characterized by double or zero oc- 
cupancy, which could be detected using parity measure- 
ments as those introduced in Refs [32, 33 1, and could be 
discerned from each other by the absence/presence of in- 
terference fringes in TOF [341 ] . Finally, the string-order 
of the HI phase may be studied using similar site-resolved 
measurements as those recently reported for the measure- 
ment of non-local parity order in Mott insulators [2(| ■ 

This work has been supported by the cluster of ex- 
cellence QUEST (Center for Quantum Engineering and 
Space-Time Research). T.V. acknowledges SCOPES 
Grant IZ73Z0-128058. 



[1 
[2 
[3 
[4 

[5 

[6 

ir 

[8] 

[9] 
[10 

[11 
[12 

[13 
[14 

[15 
[16 

[17 
[18 

[19 

[20 



M. Lewenstein et al, Adv. in Phys. 56, 243 (2007). 
C. Becker et al., New J. Phys. 12, 065025 (2010). 
G.-B. Jo et al, Phys. Rev. Lett. 108, 045305 (2012). 
S. D. Huber and E. Altman, Phys. Rev. B 82, 184502 

(2010) . 

A. Eckardt, C. Weiss, and M. Holthaus, Phys. Rev. Lett. 
95, 260404 (2005). 

A. Zenesini et al, Phys. Rev. Lett. 102, 100403 (2009). 

J. Struck et al, Science 333, 996 (2011). 

C. Chin et al, Rev. Mod. Phys. 82, 1225 (2010). 

A. J. Daley et al, Phys. Rev. Lett. 102, 040402 (2009). 

I. Titvinidze et al, New Journal of Physics 13, 035013 

(2011) . 

L. Bonnes and S. Wessel et al, Phys. Rev. Lett. 106, 
185302 (2011). 

Y.-C. Chen, K.-K. Ng, and M.-F. Yang, Phys. Rev. B 

84, 092503 (2011). 

M. J. Mark et aZ.. larXiv:1201.1008l 

E. G. Dalla Torre, E. Berg, and E. Altman, Phys. Rev. 
Lett. 97, 260401 (2006). 

T. Lahaye et al, Rep. Prog., Phys. 72, 126401 (2009). 

F. D. M. Haldane, Phys. Lett. A 93, 464 (1983); Phys. 
Rev. Lett. 50, 1153 (1983). 

E. Berg et al, Phys. Rev. B 77, 245119 (2008). 

M. den Nijs, and K. Rommelse, Phys. Rev. B 40, 4709 

(1989). 

T. Giamarchi, Quantum Physics in One Dimension, Ox- 
ford University Press (2003). 
M. Endres et al, Science 334, 200 (2011). 



5 



[21] A. A. Nersesyan, A. O. Gogolin, and F.H.L. Essler, Phys. 
Rev. Lett. 81, 910 (1998). 

[22] See supplementary material for details. 

[23] One may map boson problem by means of the Holstein- 
Primakoff transformation Sf = 1 — n;, = V2 — nbi 
to frustrated spin-1 XY chain (up to particle-hole 
symmetry breaking terms that are irrelevant [l7|). 

Hs = «E, [(Si ■ S i+ i) XY +j{Si ■ S i+2 )xY + D(S!) 2 ] , 

where (Si ■ Sj) X Y = Sf ■ SJ + Sf ■ S? , and D = U/4t and 
j characterize, respectively, the single-ion anisotropy and 
the next-nearest neighbor frustrating exchange. 

[24] Numerically we have unambiguously confirmed the Ising 
nature of SF-PSF transition by observing the correct 
scaling of Op close to transition point. 

[25] S. R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. Rev. 
B 48, 10345 (1993); U. Schollwock, Rev. Mod. Phys. 77, 



259 (2005). 

[26] It is notoriously hard to determine precisely the KT 
phase transition lines by this method alone. The precise 
location of the phase transition lines, and the determi- 
nation of the nature of the multi-critical points close to 
j = are not pursuit in the current work. 

[27] A. K. Kolezhuk, Prog. Theor. Phys. Suppl. 145, 29 
(2002); Phys. Rev. B 62, R6057 (2000). 

[28] T. Hikihara et al, J. Phys. Soc. Jpn. 69, 259 (2000). 

[29] T. Hikihara, J. Phys. Soc. Jpn. 71, 319 (2002). 

[30] H.J. Schulz, Phys. Rev. B 34, 6372 (1986). 

[31] T. Vekua, et al, Phys. Rev. B 76, 174420 (2007). 

[32] W. S. Bakr et al, Science 329, 547 (2010). 

[33] J. F. Sherson et al, Nature 467, 68 (2010). 

[34] M. Greiner et al, Nature 415, 39 (2002). 



Supplementary material to 
"Ultra-cold bosons in zig-zag optical lattices " 

S. Greschner, 1 L. Santos, 1 and T. Vckua 1 
1 Institut fur Theoretische Physik, Leibniz Universitdt Hannover, 30167 Hannover, Germany 



(N 
O 



C/2 

i 



o 
o 



> 

00 

m 
in 

o 



X 



In this supplementary material, we provide some addi- 
tional details on the bosonization description of the possi- 
ble sequential CSF-CMI-MI phase transitions for j > 1/4 
As mentioned in the main text, in the j > 1 regime we 
may consider sub-chains formed by the even and odd 
sites, and introduce the symmetric and antisymmetric 
fields 9± and cf>± . The interaction between these fields is 
given by the last two terms of Eq. (3) of the main text: 



/ H lnt = \d x 9 + sinv27r6>_— 2M cosv 27r</> + cosv 27n/>_ , (1) 

where the first term supports chirality and the second 
one favors a MI phase. 

Starting from U ^> t' , deep in the Mott phase of each 
sub-chain, then (cos s/2tt<P±) ^ 0, and thus the <p± fields 
are pinned in the Mott phase. In this case, to second 
order in A, one can integrate out in the partition func- 
tion d x 9 + from the first term of Eq. ([]]), obtaining the 
following contribution in the antisymmetric sector, 



A 2 
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J d 2 xd 2 y(d x + d y 6 + ) sin V2^0- (x) sin y/2^9- (y) 



where the average is performed in the ground state of 
the MI phase, where (d x 0+d y + ) is short ranged. Hence 
the leading contribution in the antisymmetric sector, af- 
ter carrying out operator product expansion, is a term 
~ A 2 cos \Z8tt9- . Note that a ~ —X 2 (d x 9-) 2 contribution 
decreasing the value of AT_ is obtained as well in the an- 
tisymmetric sector. The competition between cos y/8n9- 
and cos y/2~7T<f>- (obtained using mean-field decoupling of 
the last term of Eq. (JTJ) in the MI phase) is resolved with 
an Ising phase transition in the antisymmetric sector with 
increasing X/U, leading to the pinning of 0_ in the new 
ground state (^/8n0_) = n, so that (sin \12~k9- ) ^ 0, 
driving the symmetric sector into a state with finite topo- 
logical current, (d x 9 + ) ^0. 



The simplest scenario to establish the Ising phase tran- 
sition is for A'_ = 2. In that case, performing a mean- 
field decoupling of the second term in Eq. ([T]), the anti- 
symmetric sector is governed by the Hamiltonian density, 



{d x <j>. 



(d x 



—M cosv 47r</>_ + A cos v 47r6L , 



where 6L = y/KZO-, <j>- = <j>-/y/Kl, M = 
2A4(cos v / 27r0+), and A ~ A 2 . The antisymmetric sec- 
tor can be hence described by two free massive Majorana 
fcrmions, with masses A ± M. At the Ising phase tran- 
sition, A = ±M, and the mass of one of the Majorana 
fermions vanishes [lj]. 

However, the Mottness of the ground state after the 
chirality gets long-range ordered, (\/27T0_) 7^ 0, does not 
necessarily disappear immediately, due to the possibil- 
ity of a relevant contribution in the symmetric sector, 
cos \/8tt(/) + , for K + < 1 which stems after integrating 
out the cos s/2~n(f)- in the last term of Eq. (jTJ) in the state 
with pinned 9-. Note that in the CMI state (0_) ^ 0, 
(d x 9+) 7^ 0, and also (<j>+) 7^ 0. Further decreasing U, at 
K + = 1, the CMI phase (K + < 1) disappears at a KT 
phase transition in favor of the CSF phase (K + > 1). 

Thus, our bosonization analysis, for j 1, suggests 
the possibility of two consecutive phase transitions with 
increasing U starting fom the CSF phase, first a KT tran- 
sition from CSF to CMI, followed by the Ising phase tran- 
sition from CMI to MI. 
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